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Se formulo un modelo analitico que permitio estimar el tiempo de iluminacion solar sobre la Tierra para 
cualquier fecha del ano y cualquier latitud, el modelo tomo la oblicuidad de la ech'ptica constante, los rayos de 
luz paralelos, la Tierra esferica y el movimiento de la Tierra circular uniforme, tambien mostro un contexto 
de la astronomi'a para ensenar ffsica basica. Se relaciono el movimiento de la Tierra alrededor del Sol con el 
movimiento del piano de luz proyectado sobre la Tierra, luego se dedujo la zona iluminada para una latitud 
dada y se calculo cl tiempo de iluminacion mediante el movimiento circular uniforme de rotation terrestre. 
El modelo se confronto con resultados numericos de la Geoscience Australia Agency hallando como error 
porcentual maximo 1,6%, el cual se atribuyo principalmente a la discrepancia entre la trayectoria real elfptica 
y la trayectoria circular tomada en este modelo. Asf, sin hacer uso de la trigonometrfa esferica se obtuvo un 
modelo analitico que explica de manera muy aproximada el tiempo de iluminacion solar en cualquier epoca 
del ano y para cualquier latitud, el modelo brinda un contexto autentico para estudiar algunos aspectos de la 
ffsica basica. 

Palavras-chave: modelo, astronomi'a, ensenanza de la ffsica, metodos de ensefianza. 

We proposed an analytical model for the calculus of illumination time of the Earth for any time of year 
and any latitude, this model assumes the obliquity of the ecliptic as constant, the light beams as parallels, the 
Earth as spherical, the movement of translation of Earth as uniform circular, also this model showed a context 
of the astronomy whereby the teachers can teach the basic physics.lt was built through a relationship between 
the movement of translation and of rotation of the wave front light, then we found the of illumination zone on 
the Earth and the illumination time is estimated in a particular latitude with the uniform circular movement 
of Earth. Present model was confronted with the numerical results of the Geoscience Australia Agency and it 
is found a maxim perceptual error of 1,6%, this value was assigned primarily to the difference between the 
circular trajectory, in this model, and the elliptical trajectory that is the real. Without the use of spherical 
trigonometry was obtained an analytical model that estimates very close the solar illumination time at any 
time of year and any latitude on earth, the model provides an authentic context for studying basic aspects of 
physics. 
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1. Introduction 

Recientes investigaciones en ensenanza de la ffsica 
nan venido mostrando que la astronorma trae consi- 
go contextos autenticos de alta motivation en los es- 
tudiantes que han permitido desarrollar diversos con- 
tenidos fisicos a diferentes niveles de formation, tales 
contextos abarcan topicos como la Tierra, la escala del 
sistema solar, la determination de distancias a los pla- 



netas y a las estrellas, la ley de Hubble, la ley de Wien 
y las curvas de rotation en galaxias espirales, entre 
otros [1-11]. Asf mismo, estos contextos han permiti- 
do introducir contenidos fisicos tales como dinamica, 
cinematica y optica. En [2] se muestra la experiencia 
de buscar planetas extrasolares a traves de la adquisi- 
cion de datos propios mediante el control remoto de un 
telescopio disehado para tal fin; en [3] se expone un es- 
tudio que revela algunas ideas erroneas sobre la trayec- 
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toria de los planetas alrededor del sol; en [4, 5] se estu- 
dia el movimiento planetario a traves de experimentos 
virtuales con base en argumentos ffsicos elementales e 
intuitivos que permiten desarrollar discusion sobre la 
ley de gravitation, la tercera ley de Kepler y metodos 
numericos de implementation simple; en [7] se estu- 
dian curvas de rotacion en galaxias espirales ;en [8] se 
muestra como obtener la excentricidad de la orbita te- 
rrestre usando, principalmente, las leyes de Kepler y 
un instrumento de observation sencillo; en [9] se repa- 
san algunos metodos geometricos mediante los cuales 
se estima las distancias Tierra-Luna y Tierra-Sol, asf 
como los diametros del Sol y la Luna en relation al 
radio de la Tierra; en [10] se aprovecha el contexto de 
la ley de Hubble para promover la habilidad de la in- 
terpretation conceptual en graficas de velocidad contra 
distancia y en [11] se muestra como una presentation 
desde la estadi'stica de fotones de la ley de radiation 
de cuerpo negro es mas eficaz, desde el punto de vista 
pedagogico, que el tratamiento traditional que se hace 
generalmente en los libros de texto. 

Dentro del contexto anterior, este trabajo muestra 
como el fenomeno cotidiano del dfa y la noche puede 
usarse para estudiar aspectos basicos de fi'sica como 
factores de conversion, movimiento circular uniforme 
(MCU) y optica de rayos. El problema concreto de 
estudio consiste en responder la pregunta: les posi- 
ble estimar mediante un modelo analftico basado en 
conceptos ffsicos elementales el tiempo de iluminacion 
solar sobre la Tierra para cualquier fecha del aho y 
en cualquier latitud? Para responder esta pregunta se 
desarrolla un modelo analftico que deja de lado algu- 
nos hechos astrofisicos que no son relevantes durante 
el transcurso de pocos anos, asf las caracteristicas del 
modelo son la siguientes: toma en cuenta la correction 
estandar de la refraction de la luz por la atmosfera y el 
tiempo diario de iluminacion debido al movimiento de 
traslacion de la Tierra alrededor del Sol; desprecia los 
efectos tanto de precesion como de nutation; considera 
la Tierra completamente esferica; asume el movimien- 
to de traslacion como un MCU y toma los rayos de luz 
solar paralelos. 

Para la formulation del modelo se establece el siste- 
ma de referenda inercial con uno de sus ejes paralelo 
al eje de rotacion de la Tierra (considerado estatico 
respecto a las estrellas fijas) y los demas ejes tambien 
anclados a estrellas fijas. Ahora, se establece un vinculo 
entre el MCU de la Tierra alrededor del Sol y el mo- 
vimiento del frente de onda piano de luz solar el cual 
gira sobre la esfera terrestre, luego se delimita la zona 
de la Tierra iluminada en cualquier fecha del aho a una 
latitud arbitraria y se aplican conceptos de cinematica 



del MCU de rotacion de la Tierra para calcular el tiem- 
po que un punto terrestre estarfa iluminado. Mediante 
las ecuaciones obtenidas es posible conocer de mane- 
ra muy aproximada el tiempo de iluminacion para una 
fecha particular del aho en todo el mundo, conocer la 
fecha de ocurrencia del Sol de media noche para cual- 
quier latitud y mostrar diagramas de simetrfa utiles a 
la hora de interpretar hechos ffsicos mediante la lectura 
de graficos. 

El modelo permite realizar analisis grafico del tiem- 
po de iluminacion durante un aho para una latitud 
determinada, estas graficas permiten realizar una con- 
frontation sencilla del modelo con resultados numeri- 
cos obtenidos de simulaciones que toman en cuenta los 
efectos de otros fenomenos astrofisicos en particular de 
la trayectoria elfptica de la Tierra, al confrontar el mo- 
delo con los resultados arrojados en el simulador de la 
Geoscience Australia Agency se halla un error porcen- 
tual maximo en la grafica mencionada muy cercano a 
1,6% error que se atribuye fundamentalmente, y de 
acuerdo con la tercera ley de Kepler, a la diferencia 
entre las trayectorias elfptica y circular. 

Este trabajo muestra, sin hacer uso de la trigonome- 
trfa esferica, una manera de pensar ffsicamente sobre 
el hecho cotidiano de la noche y el dfa, a traves de un 
modelo analftico basado fundamentalmente en MCU 
que a traves de sus resultados permite explicar de ma- 
nera muy aproximada el tiempo de iluminacion solar 
en cualquier epoca del aho y para cualquier latitud, el 
modelo deja un escenario cautivante dentro de la as- 
tronomfa para estudiar algunos aspectos basicos de la 
ffsica como el MCU y la optica de rayos paralelos. 

Este artfculo esta organizado de manera siguiente: 
en la section (2) se deduce el tiempo de iluminacion 
terrestre de acuerdo con lo planteado anteriormente, 
en la section (3) se hacen correcciones por traslacion y 
refraction, en la section (4) se comparan los resultados 
del modelo con resultados numericos, en la section (5) 
se proponen algunas estrategias didacticas para poner 
en practica los resultados de este modelo y en la section 
(6) se muestran las conclusiones. 

2. Tiempo de iluminacion solar so- 
bre la Tierra 

Se sabe que debido a la rotacion de la Tierra sobre su 
propio eje esta no posee una geometrfa que corresponda 
por completo a una esfera sino mas bien a un elipsoide 
[12-14]. Debido a la geometrfa elipsoidal de la Tierra y 
a la interaction gravitational con la Luna, el Sol y los 
demas cuerpo celestes que componen el sistema solar 
el movimiento de traslacion y rotacion terrestre varia 
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respecto a lo que se espera del estudio de un sistema donde r Q = pq y 
Sol- Tierra completamente aislado. Sin embargo, estos 
efectos no se toman en cuenta en este trabajo debido a 
que sus efectos son despreciables durante el transcurso 
de un ano [13, 15]. 

A lo anterior se ariade que la luz solar se toma con 
la aproximacion de rayos paralelos, esto de que 

observaciones bien establecidas muestran una desvia- 
cion del paralelismo de aproximadamente 0,5°,que es 
producida por la refraccion de la luz al atravesar la 
atmosfera [16]. 

PNT 



tan e = 



h ' 

Al eliminar h de (1) y (2) se obtiene 
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— = tan^tane. 
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En la Fig. 2 se observa la Tierra aun en solsticio de 




Figura 1: Corte meridional de la Tierra paralelo a los 
rayos de luz provenientes del Sol para el solsticio de 
invierno en el hemisferio norte. 



En la Fig.l se muestra un corte meridional de la 
Tierra paralelo a los rayos de luz provenientes del Sol 
para el solsticio de invierno en el hemisferio norte, se 
observa que el angulo e formado entre el piano M del 
frente de onda de la luz solar y el eje de rotacion de 
la Tierra L es la oblicuidad de la ech'ptica, que es el 
angulo formado por el eje de rotacion de la Tierra y el 
piano de traslacion alrededor del Sol. Esto implica que 
un observador en la superficie de la Tierra a una latitud 
4> tendra una trayectoria circular con radio r = qs, una 
velocidad angular oj y un periodo de rotacion T = ^ 
que corresponden a la velocidad angular y al periodo 
de rotacion terrestre como muestra la Fig.l.Asf mismo, 
de la la Fig.l se puede obtener las ecuaciones 



tan <p = — 
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Figura 2: Corte de la Tierra a latitud (p vista desde 
PNT. 



invierno vista desde el polo norte terrestre (PNT), se 
observa que la zona iluminada es el arco de circunfe- 

rencia AB con angulo 9 y radio r = Aq. Ademas, la 
geometrfa evidenciada que el angulo 9 se puede expre- 
sar como cos (|) = lo cual permite expresar con 
ayuda de (3) el angulo 9 a traves de la ecuacion 



cos 



tan <A tane. 



(4) 



Ahora, si se supone que el piano M no rota significati- 
vamente al rededor del punto p durante una revolucion 
terrestre el tiempo de iluminacion to en el solsticio se 
puede obtener al calcular el tiempo que un observador 
sobre la superficie terrestre tarda en barrer el angulo 9 
con velocidad angular u, siendo u = j-. Asf se obtiene 



to 



T 

— cos -1 (tan < 

7T 



tan e) 



(5) 



(1) 



El tiempo tjy que el observador permanece en la zona 
no iluminada se puede calcular a partir de (5) y es 
simplemente tjv = T — to- 



3 



Al tomar como origen de tiempo el solsticio de in- 
vierno en el hemisferio norte la traslacion de la Tierra 
alrededor del Sol se ve reflejada en la rotation del piano 
M alrededor del punto p, ver Fig. 3. Si se considera el 
movimiento de traslacion como un MCU con velocidad 
angular Q entonces el piano M girara sobre p barriendo 
un angulo Qt a partir de la position initial, como mues- 
tra la Fig.3. Notese en la Fig. 3 que el triangulo qAB 



latitud la ecuacion 

4>max = tan" 1 ( ±l , (8) 

ytane cos (Ut) J 

que demarca el comienzo de Sol de media noche y la 
noche polar, esto sin tomar en cuenta el Leve incre- 
mento en el tiempo de iluminacion debido a refraction 
de la luz y traslacion de la Tierra. 




Figura 3: Position angular del piano M en tiempo t. 

es isosceles y por lo tanto al usar la suma de angulos 
internos se tiene = ir — 2a, al aplicar el teorema del 
seno al triangulo qqB se encuentra r sin a = r Q cos (Q,t) 
y usar la ecuacion (3) se encuentra 

T 

tp = — cos" 1 (tan</> tan e cos (Hi)) . (6) 

7T 

Si se considera la traslacion de la Tierra alrededor del 
Sol como un MCU con << u, se desprecia la ro- 
tation del piano M del frente de onda de la luz solar 
y tambien se desprecia la refraction de la luz por la 
atmosfera la ecuacion (6) permite estimar el tiempo 
de luz solar durante una revolution de la Tierra para 
cualquier fecha del ano y cualquier latitud. 

La ecuacion (6) tiene sentido solo si se cumple la 
condition 

|tan0 tane cos($7t)| < 1, (7) 

para que pueda evaluarse la funcion arcocoseno. Esta 
condition demarca la latitud maxima 4>max que puede 
ser evaluada en (6) de acuerdo con el tiempo t transcu- 
rrido. Asf, de acuerdo con (7) se tiene para la maxima 



3. Correccion de tp por traslacion y 
refraccion 

Durante una revolution de la Tierra sobre su eje el 
piano M no permanece en reposo, como se supuso, sino 
que rota un poco alrededor del punto p. En la Fig. 4 se 
muestra el angulo de la zona iluminada initial 6i, es 
decir justo cuando un observador sobre la superfitie 
terrestre llega al punto A y se muestra tambien un 
angulo Of de la zona iluminada cuando el mismo ob- 
servador justo alcanza el piano M por su otro extremo 
en el punto D. De acuerdo con esto el observador ba- 
rre un angulo (3 adicional, que es preciso calcular para 
determinar el tiempo adicional de iluminacion. 

En la Fig.4 se observa la relation entre los angulos 
7, at y Qt 

7 + Qi + Sir = vr, (9) 

donde el tiempo r se refiere al tiempo transcurrido des- 
de que un observador ligado a la Tierra con latitud <f> 
coincide con el piano de M en el punto A hasta que 
nuevamente coincide en el tiempo t + r, siendo t el 
tiempo transcurrido desde el solsticio de invierno en el 
hemisferio norte hasta que observador y piano coinci- 
den en el punto D. 

Igualmente en la Fig.4 se puede observar tambien 
que el angulo /3 adicional que barre el observador esta 
relacionado con los angulos 7 y «j mediante 

/3 + 7 + a/ = 7r. (10) 

Ademas, de los triangulos qAB y qCD se encuentran 
las relaciones 

2o i + i = 7T. (11) 

2a f + f = -iT. (12) 

Ahora, al combinar las ecuaciones (9), (10), (11) y (12) 
se encuentra 

P = + «r, (13) 

donde 9 t = 0(t) y f = 0(t + r) con 

0(t) = 2 cos" 1 (tan 4> tan e cos (tit)) . (14) 
Cuando la Tierra gira alrededor del Sol la zona ilumi- 
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Figura 4: Zonas iluminadas al principio y al final de 
una fraction de revolution terrestre. 



nada cambia de Oi a Of lo que implica que un obser- 
vador que coincidida con el piano de luz inicialmente 

en A barre un arco AD y emplea un tiempo r que se 
concreta justo cuando el observador alcanza el piano 
de luz en D. De manera que 



+ /3 = LOT, 



(15) 



siendo r el tiempo que el observador esta iluminado. 
Ahora, como to el tiempo que el observador esta ilu- 
minado cuando no se toma en cuenta la pequena tras- 
lacion que la Tierra sufre mientras hace una revolution 
sobre su eje At = r — Tj es el tiempo adicional de ilu- 
minacion que debe soportar el observador si se toma 
en cuenta esta pequena traslacion. 

Para calcular r se reemplaza la ecuacion (13) en (15) 
y se encuentra 

0f + 0i , 1R * 

r= 2(^y- ( 16 ) 

En (16) se espera que Of ps 0{ dado que to « u, lo 
que sugiere hacer una expansion en serie de Taylor de 
0(t) alrededor de y tomar la aproximacion con los 
primeros dos terminos de la serie, asi se obtiene 



con 



Of = 0i + 2r/T 



fitane tan0 sin(ilt) 
a/1 — tan 2 e tan 2 (j) cos 2 (tot) ' 



(17) 



(18) 



que tambien cumple la condition (7) y los terminos 
de coseno y seno junto con los ordenes de magnitud 
garantizan que (18) siempre permanece finito y menor 
que 1^. 

Para encontrar la correction por traslacion At se 
reemplaza la ecuacion (17) en (16), obteniendo 



r = 



Oi 



LO 



to 



lo que conduce a 



At 



to (to — to — rj) ' 



(19) 



(20) 



que es la correcion por traslacion buscada y resulta ser 
de w 2 m . 

De otro lado, se sabe que cuando los rayos de luz 
penetran en la atmosfera sufren refraction que aumen- 
ta el tiempo de iluminacion, asi en la Fig. 2 los rayos 
de luz paralelos a la semirecta ps que llegan tanto al 
punto A como al B convergen un poco hacia el punto 
q. De esta manera un observador comienza a ser ilu- 
minado un poco antes de llegar al punto A y sigue un 
poco iluminado despues de pasar por el punto B. En la 
literatura suele tomarse la portion de angulo que ba- 
rre un observador bajo los efectos de la luz refractada 
como Or ~ 9, 95 x I0~ 3 rad [14, 17]. De lo anterior se 
entiende que este fenomeno genera un tiempo adicional 
de iluminacion dado por que es aproximadamente 

De acuerdo con las anteriores consideraciones la 
ecuacion (6) toman la forma 



T 

t D = At + — [20 R + cos -1 (tan^ tan e cos (tot))] . 

(21) 



4. Modelo vs resultados numericos 

De los datos numericos ofrecidos por la simulation 
de la Geoscience Australia Agency puede obtenerse re- 
sultados que calculan el tiempo que un observador per- 
manece iluminado durante una revolution de la Tierra. 
Los aspectos mas relevantes para la confrontation del 
modelo analftico que toma en cuenta el modelo nume- 
rico de Geoscience Australia Agency son la trayectoria 
elfptica de la Tierra alrededor del Sol y la refraction 
estandar de la luz debida a la atmosfera [17]. 

En la Fig. (5) se muestra el tiempo de iluminacion 
to en minutos como funcion del tiempo t transcurrido 
en dfas desde el comienzo del solsticio de verano de 
2009 hasta justo el comienzo del solsticio de verano 
de 2010, se observa mayor diferencia entre los tiempos 



5 





850 




825 




800 




775 


5 


750 




725 


-w 


700 




675 




650 




625 




600 















— ^ 




nulado 










Sii 


9 



























































































61 



122 



183 
t(d) 



244 305 366 



1440 



1200 



960 



720 



480 



240 









1 1 

MCU 

Simulado • 


\ 





























































-90 -60 -30 30 60 90 



Figura 5: Modelo analftico con MCU vs simulacion pa- 
ra 4> = 30° y e = 23,5° entre el 21 de junio de 2009 y 
el 21 de junio de 2010. 



Figura 6: Modelo analftico con MCU vs modelo nume- 
rico para el solsticio de invierno en el hemisferio norte. 



obtenidos con la simulacion y el modelo analftico desde 
mediados de febrero hasta los primeros dfas de octubre 
epoca del aho en la cual la Tierra esta ceca al perihelio, 
que ocurre aproximadamente el 4 de enero [18]. Por el 
contrario se observa un buen acuerdo entre mediados 
de octubre y los primeros dfas de febrero epoca del 
ano en que la Tierra esta cerca a su afelio, que ocurre 
aproximadamente el 4 de julio [19]. La tercera ley de 
Kepler implica que cuando la Tierra se encuentra en el 
perihelio su rapidez es mayor que cuando se encuentra 
en el afelio lo que trae como consecuencia que el tiempo 
adicional de iluminacion por incremento de traslacion 
aumente ya que el piano M de la Fig.l se mueve mas 
rapido sobre el punto q y por tanto un observador sobre 
la Tierra tarda mas tiempo en alcanzar el piano de luz 
solar. Por el contrario cuando la Tierra se encuentra en 
el afelio el piano M se mueve mas lento produciendo 
que el tiempo adicional de iluminacion por traslacion 
sea menor. Consecuentemente en cecanias del perihelio 
la simulacion muestra resultados alejados del modelo 
analftico de trayectoria circular que corresponde a unos 
w 12,5 m que trae consigo un error relativo porcentual 
de » 1,6%. 

De otro lado la Fig. (6) muestra en el solsticio de 
invierno el tiempo de iluminacion para diferentes lati- 
tudes, se observa un muy buen acuerdo entre el modelo 
simulado y el modelo analftico. 

Diversas comparaciones entre el modelo analftico y 
la simulacion arrojan resultados similares. Asf, el mo- 
delo analftico da cuenta de manera cercana y con un 



mecanismo ffsico de explicacion al tiempo de ilumina- 
cion solar en cualquier tiempo del aho y en cualquier 
lugar de la Tierra. 

5. Algunas sugerencias didacticas 

De acuerdo con los resultados obtenidos la ecuacion 
(21) permite estimar el tiempo de iluminacion para una 
fecha particular del aho y para un lugar particular en 
la Tierra caracterizado por su latitud, tomando correc- 
ciones por refraccion y traslacion. Con este resultado 
es posible atacar problemas de la vida cotidiana tales 
como la estimation del tiempo de iluminacion solar en 
un pais particular; apoyar una explicacion ffsica cuan- 
titativa de las estaciones del aho; ofrecer un espacio 
para conceptualizar el MCU ejemplihcando a traves 
de los movimientos de traslacion y rotation de la Tie- 
rra; usar esquemas como el mostrado en la Fig. 3 para 
plantear una discusion sobre la refraccion de la luz por 
la atmosfera; discutir los fenomenos del Sol de media 
noche y la noche polar, en tre otros. Asf mismo, se pue- 
den plantear situaciones en las cuales dada una fecha 
particular del aho se pregunta por la duration de la 
iluminacion en diferentes pafses y ofrecer explicacion 
ffsica de las diferencias presentadas. Resulta particu- 
larmente interesante usar el modelo para estimar la 
fecha de initio de los dfas y noches polares, para es- 
to es necesario realizar analisis grafico de la ecuacion 
(8) que muestra la maxima latitud (p ma x que es posible 
evaluar en las ecuaciones (5) y (21) como funcion del 
tiempo y que a su vez permite estudiar las zonas de 
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Figura 7: Latitud para el Sol de media noche y la no- 
che polar con el transcurso del aho sin correciones por 
refraccion ni traslacion. 



la Tierra donde comienza el sol de media noche o la 
noche polar. 

En la Fig. 7 se muestra la evolution de la latitud ma- 
xima con el transcurso del aho para los hemisferios 
norte y sur, asi por arriba de la lmea de puntos oscu- 
ros habra latitudes que estan en noche polar para una 
fecha determinada del aho y por debajo de la lmea de 
puntos claros habra latitudes en las cuales habra dfa 
polar. 

Ademas, el modelo puede usarse para estudiar al- 
gunos aspectos referidos a la iluminacion de planetas 
tonto en el sistema solar como en otros sistemas sola- 
res. 

6. Conclusiones 

Se demostro que es posible formular un modelo ana- 
h'tico con base en la hipotesis de un MCU de la Tierra 
alrededor del Sol que permite estimar el tiempo de ilu- 
minacion para un tiempo cualquiera despues del sols- 
ticio de invierno en el hemisferio norte a cualquier lati- 
tud, que concuerda bien con el modelo numerico de la 
Geoscience Australia Agency y se realizaron correccio- 
nes por traslacion y refraccion siendo respectivamente 
de ~ 2 m y de ~ 4 m 5 s , el modelo brinda un escenario 
rico en fenomenos terrestres relacionados con la astro- 
nomfa que permiten estudiar algunos aspectos de la 
fi'sica basica. Queda como problema sugerido aplicar 
el modelo a planetas extrasolares con diversas caracte- 
n'sticas para estudiar diferentes aspectos relacionados 



con el tiempo de iluminacion sobre sus superficies. 

Agradecimientos 

Los autores agradecen a los profesores Fabian Ga- 
lindo, Giovanni Cardona y Guillermo Avendaho por la 
lectura y las sugerencias realizadas sobre el manuscrito; 
al Departamento de Ciencias Naturales de la Univer- 
sidad Central por el apoyo y el tiempo asignado a la 
investigation y al CICATA del IPN de Mexico por su 
continua colaboracion. 

Referencias 

[1] C. H. Holbrow and P. Shaffer, Am. J. Phys. 80, 
362 (2012). 

[2] R. Goulda and S. Sunbury, Am. J. Phys. 80, 445 
(2012). 

[3] F. Santana y P. Talero, Lat. Am. J. Phys. Educ. 
5, 492 (2011). 

[4] F. Santana y P. Talero, Lat. Am. J. Phys. Educ. 
5, 508 (2011). 

[5] P. Talero Experimentos Virtuales de Fisica: Me- 
cdnica (Editorial Universidad Central, Bogota, 
D.C., 2009). 

[6] F. Santana Estudio de la mecdnica celeste del viaje 
a Mate con experimentos virtuales (Universidad 
Distrital Francisco Jose De Caldas, Bogota, D.C., 
2011). 

[7] O. D. Pavioni y F M. Ortega, Lat. Am. J. Phys. 
Educ. 3, 320 (2007). 

[8] O. D. Pavioni y F M. Ortega, Lat. Am. J. Phys. 
Educ. 3, 320 (2007). 

[9] P. Nunez, S. E. Calderon y S. Gil, Lat. Am. J. 
Phys. Educ. 3, 329 (2009). 

[10] C. S. Wallacea and E. Pratherb, Am. J. Phys. 80, 
324 (2009). 

[11] J. M. Marra and F. P. Wilkinb, Am. J. Phys. 80, 
399 (2009). 

[12] H. Karttunen,P.Kroger,H.Oja, M.Poutanen and 
K. J. Donner Fundamental Astronomy (Springer, 
New York, 2003). 

[13] R. M. Green Spherical astronomy (Cambridge 
University press, London, 1985). 



7 



[14] J. G. Portilla Elementos de astronomia de posi- 
tion (Editorial Universidad Nacional de Colom- 
bia, Bogota, D.C., 2001). 

[15] M. Alonso y E. J. Finn Fisica Volumen I: Mecdni- 
ca (ADDISON- WESLEY IBERO AMERICANA, 
U.S.A, 1976). 

[16] P. G. Hewitt Fisica conceptual (PEARSON 
ADDISON- WESLEY, Mexico, 2007). 

[17] http:/ /www. ga.gov. au/geodesy/astro/sunrise.jsp# 
intzone, consultado el 07/17/2012. 

[18] http : / / ciencia. nasa. gov / science-at - 

nasa/2001/ast04jan_l/, consultado el 
07/23/2012. 

[19] htt p : / / ciencia. nasa. gov / science- at - 

nasa/2002/02jul_aphelion/, consultado el 
07/23/2012. 



8 



